Minimum Electrical and Thermal Conductivity of Graphene: 
A Quasiclassical Approach 
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We investigate the minimum conductivity of graphene within a quasiclassical approach taking into 
account electron-hole coherence effects which stem from the chiral nature of low energy excitations. 
Relying on an analytical solution of the kinetic equation in the electron-hole coherent and incoherent 
cases we study both the electrical and thermal conductivity whose relation fuUfills Wiedemann-Franz 
law. We found that the most of the previous findings based on the Boltzmann equation are restricted 
to only high mobility samples where electron-hole coherence effects are not sufficient. 



Introduction. Single graphite layers (graphene) have 
been found in the free state only recently [Ij], and their 
transport properties have immediately attracted much 
attention from both experimental [1, [3^ 1, S S Hi 
and theoretical [1, 0, [ul, [H, IH M, O, Hi investiga- 
tors. The reason of such an explosive interest is a number 
of very unusual transport properties including (i) a non- 
vanishing electrical conductivity even at zero carrier con- 
centration (minimum conductivity phenomena), and (ii) 
independence of this minimum conductivity of tempera- 
ture. Besides these unconventional transport properties 
many related phenomena have been studied in graphene 
such as weak localisation [l5|, the Klein paradox fid\ 
etc. The remarkable electronic properties of graphene 
are usually attributed to the particular spectrum of ex- 
citations [17] which consists of two conical bands and is 
described by a two-dimensional analog of the relativistic 
Dirac equation. For a review concerning the history, fab- 
rication, fundamental properties, and future applications 
of graphene we refer to the recent article [l^ . 

To investigate the transport properties (i) and (ii) of 
graphene several different approaches have been applied 
including the Kubo formalism [ol, direct calcula- 

tions of the transmission probability for ballistic samples 
11, 121, and Boltzmann equation [3, The latter ap- 
proach looks at the first sight inapplicable for investiga- 
tion of the minimum conductivity since the quasiclassi- 
cal description is expected to fail at low Fermi energies 
Ep as soon as EpT{Ep) becomes comparable with h. 
The carrier momentum relaxation time r, on the other 
hand, diverges at Ep = Q for short-range scatterers 
studied here as well. Thus, the product EpT[Ep) de- 
pends at Ep ^ on the scattering parameters rather 
than the carrier concentration and can acquire values 
much lager than h even at zero doping. In this Letter 
we solve the kinetic equation for Dirac fermions includ- 
ing off-diagonal elements of the distribution function in 
the helicity basis which are strongly connected with the 
electron-hole coherence or Zitterbewegung effects stud- 
ied recently in graphene by Auslender and Katsnelson 
[l^ . We find that the off-diagonal elements essentially 
contribute to the minimum conductivity in low mobil- 



ity samples, and, thus, the conclusions obtained in [l3| 
are restricted to only quite perfect graphene sheets. To 
discuss our findings we adduce both the electron-hole co- 
herent and conventional solutions of the kinetic equation. 
As an application the electrical and thermal conductivity 
is calculated. 

Preliminaries. The carriers in the 7r-system of 
graphene near half filling can be described by the Dirac 
Hamiltonian 



H = hVQ^axkx + CFyky), 



(1) 



where vq « lO^ms"^ is the effective "speed of light", 
ax,y are the Pauli matrices, and k is the two-component 
particle momentum. The eigenstatcs of ([1} have the form 



(2) 



where tanO = ky/k^, and the energy spectrum reads 
Ek± = zthvok. The velocity matrix in the basis ^ is 



V 

Vq 



COS 9 — isin^ 
i sin — cos 6 



sin0 
~i cost 



I cose 
— sin i 



(3) 

Let us first consider the Boltzmann equation for charge 
carriers in the presence of impurity scatterers described 
by the effective potential 



V{v) 



qeZ 



-r/R 



(4) 



where eZ, q are the impurity atom and carrier electro- 
static charge respectively, and R is the screening radius. 
The potential ([4]) differs from its conventional short-range 
(5-function approximation Q by the additional fitting pa- 
rameter R. As we shall see below, the minimum conduc- 
tivity value is governed by both the impurity concentra- 
tion and screening radius whereas the carrier mobility 
turns out to be i?-independent. It seems to be necessary 
to introduce such a parameter in order to explain the ex- 
perimental picture [8| where two samples made out of the 
same graphene flake (i. e. having equal mobility) demon- 
strate essentially different minimum conductivity values. 
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This difference is attributed to the screening parameter 
R in our model. 

Electron-hole incoherent solution. In hnear order in 
the homogeneous electric field E the Boltzmann equation 
reads 



jr \ coll 

dfK\ „ 



(5) 



where we have divided the distribution function (k) — 
/°(-E'fcK)+/^(k) into an equilibrium contribution /'^{Ekn) 
and a nonequilibrium part /^(k), and (with k € {±}) 
are the diagonal elements of the velocity operator in the 
helicity basis, cf. Eq. ([3]). Assuming elastic scattering ful- 
filling the microreversibility condition, the collision term 
can be written as 

(f')'°" = E/ ^Mk.,kV)[/.,(k')-/K(k)]}, 

(6) 

where the scattering probability z/;(kK, k'/t') can be easily 
found from Fermi's golden rule, 



7 0{tjkK 



KK COS 



Ek'K') X 
'0) 



(7) 



1 + i?2 + _ 2kk/ COs{e' - 0)] 



u;(kK, k'^t') = 



Here Vq ~ TrqeZyN is an effective scattering poten- 
tial with N being the impurity concentration, and the 
fourfold (valley and spin) degeneracy is taken into ac- 
count by the factor 4 in the integrand. The exact ana- 
lytical solution of Eq. (O can be given in terms of the 
nonequilibrium part of the distribution function = 
gEv,T(fc) [-dfiEk^)/dEk^] with T{k) given by 



-{k) 



Wo 



1 

k'mw^ 1 . 
1 h^vo 



41,4 



Rk < 1. 



(8) 
(9) 



Electron-hole coherent solution. So far we have ne- 
glected the off-diagonal elements of the distribution func- 
tion, an valid approximation if the decoherence of the 
single-carrier state is a fast process compared to its re- 
laxation described by the above equations. In general, a 
particle described by the Hamiltonian ([T]) can not only 
be in one of the states ^'k+ or 'I'k- but in an arbitrary 
superposition of them. Therefore, generalising the above 
considerations, the distribution function is 2 x 2 nondi- 
agonal matrix /(k), and the kinetic equation contains 

the commutator jr H, /(k) , which drops out if only the 

diagonal elements of /(k) (with respect to the helicity 
basis) are retained. In the linear response regime the 



kinetic equation explicitly reads 




_ /i2 {Ek+ — E^] 

h\h2{Ek+-E^) 



-qE 



-vii 

V2 



■IE. 



1 { fO _ fO 



Vl2 ( fO _ fO ^ 
2Ek+ \JEk- JEk+) 

df(E^-V 



-V22 



dEk 



(10) 



with the collision term given by the generalised expres- 
sion |2C 




coll 



cPk' 



k' K,' 



EkK.) 



+S{Ek^^'^ - i?fe«)]X«f;, (fc') - SiEk^, 

x[i^:f:i/.'K,(fc) + <:;/<«i(A;)]}, 

being 



Ek'K{) 

(11) 



1 



kk'kik[ 



e'~e) 



1 + i?2 [k^ + k'^ - 2kk' cos{e' - 0)] 



(12) 



Note, that the collision term includes off-diagonal ele- 
ments of / and in that way significantly enhances the 
complexity of the kinetic equation. However, by some- 
what more tedious calculations one can again construct 



an analytical solution to Eq. 
terms given by 



TU)l with nonequilibrium 



1 

'2^ 



/l2 



gEviir(fc) 

' dfjEk-) 

dEk- 
gEvi2T(fc) (i 



1 + — 



df{Ek, 



2aEi 



k+ 



dEk+ 



(13) 



2k) r 1 



1 + 2iEk+T{k)/h { Ek+ 

df{Ek+) df{Ek-)- 



fp -fp 



dE, 



k+ 



dEk- 



(14) 



and /22, /21 can be obtained from Eqs. p^[H)) just ex- 
changing the indices belong to Ek and v accordingly. 
Here we have introduced the novel electron-hole inco- 
herence parameter a = 4i?^_|^T^(fc)/fi,^. In order to sim- 
plify the solution we use T(fc) given by Eq. thus, 
a = H'^Vq/AR'^Vq is independent on k. In the limit case 
of weak scattering (a ^ 1) the diagonal elements fl^ 
and /22 are the same as in the electron-hole incoherent 
case and given by k € {±}, whereas the off-diagonal 
elements /j^2 ^^'^ /21 '"^al and, as we shall see from 
Eq. (fT5|) . do not contribute to the current. 

Electrical conductivity. The electrical current reads 



(fk 
(27r)2 



(vii/11 + V22/22 + 2QV123/12) , (15) 
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where the factor 4 is due to the fourfold degeneracy. 

Let us first concentrate on the electron-hole incoherent 
case when a 1. Utilising Eq. ^ we find from Eq. p5|) 
for the conductivity 



CTO 



(16) 



which, in particular, does neither depend on temperature 
nor on the carrier concentration. Thus, the most strik- 
ing features of the experimental findings are reproduced: 
(i) the conductivity is not zero even at zero carrier con- 
centration, (ii) this minimum conductivity does not de- 
pend on temperature. Note, that the above result is not 
universal, i.e. the minimum conductivity at zero Fermi 
energy can change from sample to sarnple in accordance 
with recent experimental reports 0, [1, [lH. Moreover 
y/a > 4 (i.e. 1/a < 0.0625) for the vast majority of 
samples 0, Q. 

Let us take into account higher-order terms in T{Rk). 
Then instead of Eq. (|16p we have at zero temperature 



— T{kF)EF, 



(17) 
(18) 



where the carrier (electron) concentration is given by 
n — kp/n with kp — Ep/hvo being the Fermi wave 
vector. Thus, the low-temperature conductivity at low 
doping increases linearly with carrier concentration, in 
accordance with the experiments. Deviations from lin- 
ear dependency can be described taking into account k- 
dependence of the relaxation time given by Eq. ([H]). We 
emphasise, that in contrast with [8| we deal with two 
fitting parameters R and Vq which can be deduced us- 
ing Eq. (fT8| and the experimental data 0, 0, Hi- In- 
deed, from Eq. (fT8|) we can define the electron mobility as 
/i = qhvQ / (2Vq) which does not depend on the screening 
radius R. The effective scattering potential then reads 
Vb — vo\/ qh/(2ii), and for the most common samples 
with /i ranged from ^ lO'^ to 2 • lO'' cm^/(Vs) we have Vq 
covering the range from 0.1 meV to 0.06 eV. The screen- 
ing radius can be estimated from Eq. p6p assuming that 
ctq is of the order of Aq^/h. Then we have R of the order 
of hvo/Vo covering the range from 10^'^ cm (high mobil- 
ity samples) to 10~^cm (low mobility samples). Now 
one can see that the linear approximation psp in terms 



of 27ri?2„ at n - 10^^ 

cm holds only in relatively low 
mobility samples, in accordance with the experimental 
report 0, @]- The difference between cr(n) and its linear 
approximation can be seen in Fig. [T] 

To consider the low mobility samples properly, we 
should take into account the terms proportional to 1 /a in 
the solution (fT3UT4|) . Note, above all, that the quasiclas- 
sical approach is doubtful at a ^ 1 since Eprikp) H, 
and the carrier mean free path becomes comparable with 
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FIG. 1: Conductivity dependence on the carrier concentra- 
tion, cf. Ref. @]. Inset: the minimum conductivity vs. inverse 
mobility, cf. Ref. [3]. In order to preserve the validity of the 
quasiclassical description the parameters are chosen so that 
Of > 1. The electron-hole coherence correction is nevertheless 
clearly seen. 



its de Broglie wavelength. Nevertheless, an asymptotic 
dependence of the solution close to a ~ 1 can give us a 
clue to what happens in this regime. The direct integra- 
tion of Eq. (fT5|) leads to the logarithmic divergence due to 
the terms proportional to fp^_^ — fp^_ in Eqs. (fT3UT4|) . 
This problem has been solved in Ref. [3] introducing 
some ultraviolet cut off energy Ec which allows to get a 
finite value for the conductivity but still has an unclear 
physical meaning itself. In our opinion the divergence of 
the integral in Eq. (fT5|) is a clear manifestation of obvi- 
ous limitations inherent in the quasiclassical approach at 
EpT{kp) ~ H. Thus, we expect the ultraviolet cut off Ec 
to have a quantum mechanical origin. At Ep — and 
zero temperature the quasiparticle energy will fluctuate 
around E — with a variance AE related to the relax- 
ation time via the uncertainty relation AEt ^ h. Obvi- 
ously, Ec is the maximum energy uncertainty consistent 
with two subsequent scattering events ( "measurements" ) 
separated by a time interval T{kp), i. e. Ec ^ hT^^{kp). 
Then Eq. (llSp can be integrated easily, and the minimum 
conductivity for low mobility samples reads 





2 




4 




n 












a 




a 





(19) 



The difference between CTmin and ctq given by Eq. (|16p is 
shown in Fig. [T] (inset). We emphasise that Eq. (|19p could 
not be well mathematically grounded in the framework 
of our quasiclassical model because of quantum effects 
which obviously contribute to the conductivity minimum 
at a ~ 1. What is certainly true is that the additional 
term in the conductivity minimum stemming from the 
electron-hole coherence increases in low mobility sam- 
ples and partly compensates the diminution of the lead- 
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ing term given by Eq. (|16p . This mechanism might be 
responsible for the non-monotonic dependence of the con- 
ductivity minimum on the impurity concentration (in- 
verse mobihty) observed recently 

We have also studied the conductivity choosing poten- 
tial profiles different from one given by Eq. To give 
an example for hard-wall potential when ^(r) = Uq at 
r < ro, and ^(r) = at r > rg, we obtain for the con- 
ductivity the same formulas as before besides the substi- 
tutions Vq tttq^/NUo and R — > ro/2. The mobility 
becomes rg-dependent in this case that makes it difficult 
to fit our model to the experimental data Q . Most inter- 
esting is, however, (5-function shaped scattering potential 
which can be used as a model for neutral impurities or 
lattice imperfections. In this case the Boltzmann equa- 
tion ^ can be solved easily, and the conductivity does 
not depend on the carrier concentration at all, what, in 
turn, contradicts the measurements. Thus, our short- 
range potential choice given by Eq. (g]) fits best to nowa- 
days experimental data. 

Thermal conductivity. Following our method it is pos- 
sible to show that the thermal conductivity of graphene 
also has a minimal value which does not depend on the 
concentration, but depends on the temperature accord- 
ing to the Wiedemann-Franz law as was pointed out in 
Refs. 
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23l | . Indeed, in presence of the temperature gra- 
dient Vr we have at a ^ 1 in the linear response regime 
fl = {\IT/T)^^T{k){Ek. - EF)[-df{Ek.)/dEk.]. 
The further calculations of the thermal flow are very sim- 
ilar to that for the electrical current. In particular, the 
minimum thermal conductivity takes the form 



[Ef = 0) ^ C = 



T. 



(20) 



The influence of electron-hole coherence on 



can be 

described by Eq. ([19]) in full analogy with the electrical 
conductivity. 

Conclusions. We have solved the quasiclassical kinetic 
equation for carriers in a single graphene sheet including 
the off-diagonal elements of the distribution function in 
the helicity basis. The analytical solution allows us to in- 
vestigate the influence of the electron-hole coherence on 
the minimum conductivity phenomena as well as to dis- 
cover the limitations of previous studies based on Boltz- 
mann equation. We have introduced a special parameter 
a distinguishing the electron-hole coherent and incoher- 
ent regimes. It is noteworthy that a can be deduced 
directly from the minimum conductivity measurements 
since it is incorporated into (Tq in a simple way given by 
Eqs. (Uni), HH). Moreover, our approach successfully de- 
scribes the linear dependence of the conductivity above 
its minimum which is usual for low mobility samples. 
Finally, we predicted the existence of the thermal con- 
ductivity minimum which was not observed so far. 
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